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Abstract. Classical transport in a doubly connected polygonal billiard, i.e. the 
annulus square billiard, is considered. Dynamical properties of the billiard flow 
with a fixed initial direction arc analyzed by means of the moments of arbitrary 
order of the number of revolutions around the inner square, accumulated by the 
particles during the evolution. An "anomalous" diffusion is found: the moment 
of order q exhibits an algebraic growth in time with an exponent different from 
q/2, like in the normal case. Transport features are related to spectral properties 
of the system, which are reconstructed by Fourier transforming time correlation 
functions. An analytic estimate for the growth exponent of integer order moments 
is derived as a function of the scaling index at zero frequency of the spectral 
measure, associated to the angle spanned by the particles. The n-th order 
moment is expressed in terms of a multiple-time correlation function, depending 
onn-1 time intervals, which is shown to be linked to higher order density spectra 
(polyspectra), by a generalization of the Wiener-Khinchin Theorem. Analytic 
results are confirmed by numerical simulations. 



PACS numbers: 05.45.-a, 05.60.-k, 05.20.-y 
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1. Introduction. 

A fundamental problem in statistical mechanics is to understand how the reversible 
microscopic dynamics of particles may produce macroscopic transport phenomena, 
which are described by irreversible laws such as diffusion equation or Fourier law. 
In particular, the last few years have witnessed a large debate whether chaos at a 
microscopic level is a necessary ingredient to generate realistic macroscopic behaviour 
[2 [2]. In this paper we will be concerned with transport (diffusion) properties; for 
this purpose, billiards represent a class of dynamical systems ideally suited, as they 
allow both theoretical considerations and extensive numerical simulations: transport 
is typically studied by lifting the billiard table on the plane, like in the case of the 
periodic Lorentz gas [TJ [3j |H [6j. We also remark their physical significance as 
simplified models to study energy or mass transport in realistic systems, such as 
fluids, nanodevices, electromagnetic cavities, optical fibers and low density particles 
in porous media [7| 151 151 HP) . 

While the origin of normal diffusion (i.e. when the mean square displacement of 
the particles grows asymptotically linearly in time) is well understood in fully chaotic 
billiards [3 [11], in recent years, great interest has focused on the study of transport 
in polygonal billiards, which are characterized by the absence of hyperbolicity and 
dynamical chaos, in the sense of exponential divergence of nearby initial trajectories. 
Indeed in polygonal billiards, all the Lyapunov exponents and Kolmogorov-Sinai 
entropy vanish and the dispersion of initially nearby orbits is polynomial. Nevertheless 
polygonal billiards may give rise to a wide range of transport regimes, extending from 
normal to "anomalous" diffusion, in which the r.m.s. displays a non linear dependence 
on time [TSJ [El El [HI [16] . Some necessary conditions for the occurrence of normal 
transport in periodical billiard chains have been singled out: vertex angles irrationally 
related to tt, absence of parallel scatterers, existence of an upper bound for the free 
path length between collisions p"7] . 

Particle dynamics in billiards can be described in terms of an invcrtiblc flow in 
continuous time or, equivalently, of an invertible map connecting two collision points, 
which corresponds to a unitary evolution ruled by the Koopman operator. The main 
features of the dynamical transport are related to ergodic properties of the system, 
which can be formulated in terms of properties of the spectrum of the Koopman 
operator. Ergodic properties of polygonal billiards have been extensively tested by 
the analysis of the decay of correlations, such as mixing in the triangle with irrational 
angles |18j and weak-mixing as the maximal ergodic property in the right-triangle 
with irrational acute angles [TH]- The billiard flow in a typical polygon is ergodic, 
nevertheless the subset of rational polygonal billiards, in which all the angles between 
the sides are rational multiples of 7r, possess weaker ergodic properties. Some exact 
results are known for this subset: the billiard flow is not ergodic, because of a finite 
number of possible directions for a given initial condition, and it can be decomposed 
into directional flows, which are ergodic but not mixing, for a generic choice of the 
direction and for almost all initial conditions [2PJ [21] . 

The model considered in this paper is a doubly connected rational billiards, called 
the annulus square billiard, which was also studied in a quantum dynamical context 
[221 [23]. The billiard table, shown in figure [1] is formed by the plane region included in 
between two concentric squares with sides of different length. Arithmetical properties 
of the ratios between the sides of the two squares and between the components of 
the velocity vector of the particle determine different dynamical and spectral features, 
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which allow a classification of the system into different subclasses [23], reviewed in 
section |2"T21 We inspect the angle (and its absolute value), accumulated by a particle 
revolving around the hole, and provide a complete characterization of the transport 
process by the analysis of moments of arbitrary order of this observable. The statistical 
approach is supported by a spectral analysis of the system: the interdependence 
between dynamical and spectral properties leads to a theoretical prediction for the 
exponents of moments of integer order. We restrict to the class of billiards in which 
the ratios of the sides and velocity components are both irrational; in this case, for 
typical values of the parameters, the billiard flow, in a fixed direction, is weakly- 
mixing and not mixing, entailing the presence of a singular continuous component in 
the spectrum, namely supported on a set of zero Lebesgue measure [24] . By assuming 
that the singularities of the spectral measure are of Holder type [3SJ [25] , we derive an 
analytical relation connecting the exponents of the algebraic growth in time of integer 
moments with the scaling index at zero frequency of the spectral measure, associated 
to the angle. The formula is exact for even-order moments, while it is an upper bound 
for odd-order moments. This result generalizes the one obtained in [24| for the mean 
square displacement (i.e. second order moment). Furthermore, analysis of higher 
order statistics demands the introduction, in the context of dynamical system, of 
basic concepts as multiple-time correlation functions and polyspectra, which are rather 
familiar in signal analysis [27]. The theoretical relation for integer order moments is 
tested by numerical simulations, which moreover suggest that a similar estimate holds 
for moments of arbitrary order: the g-th order moment grows in time with an exponent 
vq; v is a constant function of the scaling index of the spectral measure in 0, which 
ranges between v — 1/2 ("normal" diffusion) and v — 1 ("ballistic" transport). 

A more detailed characterization of the anomalous diffusion process has recently 
been considered in many papers about diffusion in periodical polygonal billiard 
channel, in which the polygonal scatterers form a 1-dimensional periodical array 
[121 031 O US EH]. These studies constitute a further motivation to the analysis 
of transport in the square annulus billiard; indeed, to examine the dynamics of 
particles winding around the inner square is equivalent to consider the transport 
in a "generalized Lorentz gas" , i.e. an extended system, obtained by periodically 
repeating the elementary cell in both directions. In periodic channels, two different 
behaviours are identified: "weak" anomalous diffusion, when the exponent of the 
algebraic growth of moment of order q is a linear function of the order, namely vq 
with v — const. ^ 1/2 for Vg, and "strong" anomalous diffusion, when the exponent 
is a non trivial function of the order, i.e. v(q)q [Ml ISO]- In some polygonal chains, a 
transition from normal to anomalous diffusion has been found in cases with an infinite 
horizon, i.e. when the particle may travel without colliding with the walls, or, when 
the particle can propagate arbitrarily far by reflecting off only on parallels scatterers 
[15l [TBI ■ A va l ue Icr exists that discriminates between the two regimes with different 
slopes of the exponent versus q. This phase transition physically corresponds to a 
change in the balance between the ballistic and diffusive trajectories in the ensemble 
averages and infers the absence of a unique scaling law for the probability density 
of particles, which does not relax at long times into a self-similar profile; this kind 
of processes are therefore classified as " weak" self-similar processes. In some other 
cases, such as in the "zigzag model" [T71[2S], the phase transition is not present; this 
seems to be related to the geometry of the polygonal channel and, in particular, to 
the number-theoretical properties of the angles. 

The paper in organized as follows. The billiard model is described in section [2] 
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Figure 1. (Colour online) Billiard table for ii/?2 = 7r/2. A portion of a trajectory 
of a particle with initial condition (s, ip), evolving under with t > 0, is shown; 
(s, —tp) is the "time reversed" initial condition. 



and its basic spectral properties are reviewed. In section [3] we introduce the phase 
observable which identifies the diffusive process and we developed the spectral analysis 
of the signal produced by time evolution. We then consider arbitrary order moments 
of the observable; the crucial issue is to get informations on the scaling function, which 
determines the algebraic growth of the moments for long times. From the analytic 
point of view we restrict to integer order moments: higher order moments are expressed 
in terms of multiple-time correlation functions and polyspectra, which are introduced 
in section [U The analytic relation between spectral and dynamical indexes is derived 
in section [5] and numerically confirmed in section [6] 

2. The model. 

2.1. Description of the billiard table 

We consider the dynamics of a classical point particle in a doubly connected rational 
polygonal billiard. The billiard table, shown in figure [TJ is delimited by two concentric 
squares with parallel sides of length l\ and I2 (h > h)- The particle moves freely 
inside the billiard, colliding elastically with the boundary; the (preserved) modulus of 
the velocity is taken equal to 1. The accessible phase space is 3-dimensional and the 
representative point is (x,y,9): (x,y) are the Cartesian coordinates of the particles 
and 8 is the angle of the velocity, measured counterclockwise from the positive x-axis. 
The billiard flow, denoted by T*, preserves the Lebesgue measure dxdydO/ (2ir(lf — l^)). 

We may also consider the discrete time dynamics, induced by successive collisions. 
Each collision point can be parametrized by a coordinate s along the boundary and 
an angle cp, between the outgoing direction and the inner normal to the boundary; 
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we take tp positive if the outgoing velocity is given by a counterclockwise rotation to 
the inner normal. By denoting L = 2(Zi + l 2 ), the arclength s E [— L,L) takes the 
value s = —L in the left corner of side 1 and increases by moving counterclockwise 
along the sides of the billiard, labeled by increasing numbers (see figure [l}. The 
corresponding phase space is called M. = {(s,(p);s G [—L,L),(p e [— ir/2, tt/2)} and 
the Birkhoff-Poincare map is denoted by B l (t £ Z); the map B f preserves the measure 
dQ(z) = costpdtp ds /(^L) on M. 

Dynamics on such billiards is strongly influenced by number-theoretical 
properties of I2/I1 and of tani/?; in particular we may distinguish three classes |24j : 
(i) class 1: h/hjt&nip £ Q; (ii) class 2: h/h G Q, t&np e l-Q; (Hi) class 3: 
h/h G K — Q, tan ip € R — Q. In this paper we restrict to dass 3. 



2.2. Review on basic spectral properties 

As all the internal angles are rational multiples of tt, it is fairly easy to realize that 
the flow on the global phase space, or the mapping B* on M., can never be ergodic. 
However T* can be decomposed into the one-parameter family of directional flows at 
fixed 8, Tg, whose dynamics is not trivial: in particular, Tg is ergodic for almost all 
E [0,7r/2) and never mixing [5D1HT]. It is believed that the weak-mixing property 
should be enjoyed by a "generic" subset of directional flows, and numerical evidences 
for class 3 are provided in [24] . 

Analogously to the billiard flow, the map B* can be decomposed into a one- 
parameter family of components B^ along the fixed direction <p (with < <p < tt/2). 
Once the initial outgoing angle is set, only the angles ±(p, ±(-| — ip) can be met along 
the trajectory, so the phase space for any fixed foliation, i.e. M. v , is a fourfold replica 
of [-L, L): its points are denoted by z — (s, with ipj — ip — j7r/2 and j = 0, 1. 

The phase average of an observable / : M v — > C can be written explicitly as: 

f d£l(z)f(z) = ±r [ L ds ^COS^ {/( S ,^) + /( S ,-^)}. (1) 
JM^ XL J-L J=0 

Spectral properties of the system can be formulated in terms of the asymptotic 
behaviour of the correlation functions. 

Ergodicity of B^ implies that the time-averaged correlation function of a generic 
(real) phase function / : M v —> M: 



Cf" e (M)= Km ^ E f(K +tz )f(< z )> ( 2 ) 



and the phase-averaged correlation function: 



Cf (t) = (fW'f) = j dil(z) f(Blz)f(z), (3) 



coincide for i7-almost every z € M v . The unitary operator U in @ is the Koopman 
operator, associated to the map 

As directional flows or mappings are not mixing, correlation functions (of zero 
mean observables) do not decay to zero in the asymptotic limit; while the (conjectured) 
weak-mixing property guarantees that integrated correlations: 

t-i 



1 s=0 



2 
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Figure 2. (Colour online) (a) Phase averaged (upper line) and integrated (lower) 
correlation functions of the angle £, spanned by the radius vector, as a function 
of time t (measured in numbers of collisions). The billiard table belongs to class 
3 with parameters: h/h = 7r/2 and tanip = 7r/4. The dotted line, which fits 
C' nt (t), has a slope equal to — D%. (b) Estimates of generalized dimensions D\ 
(empty symbols) and D2 (full symbols) of the spectral measure, associated to 
£, for the same parameters. The dimension estimates are given by the slopes of 
straight lines: Di = 0.69 ± 0.01 and D 2 = 0.42 ± 0.02. Details are explained in 
section [6] 



vanish as t — > 00. Phase- averaged correlation function and integrated correlation 
function of the angle £, spanned by the radius vector of a particle between two 
successive collision points, are shown in figure [2fa) for a generic billiard belonging 
to class 3: while C ph {t) -fr (upper line), Cf\t) decays to zero polynomially (lower 
line) . 

By using the spectral decomposition of the Koopman operator U, we may rewrite 
© as 

Cf (t) = ± £^ f (6)e m , (5) 

which provides a direct link between the autocorrelation function of an observable / 
and the associated spectral measure. If (in the complement of constant functions) the 
spectral measure is absolutely continuous, the system is mixing. On the other side, 
weak mixing is equivalent to an empty point spectrum, apart from the eigenvalue 1. 
Therefore, the weak mixing property, without the stronger mixing property, entails 
the presence of a singular continuous component of the spectrum of the Koopman 
operator. 

Owing to the presence of a non empty point spectra, weakly mixing property is 
ruled out in the almost integrable cases, namely for Z1/Z2 S Q (class 1 and class 2); 
we will not consider such cases in the present work. 

In [24], the occurrence of a singular continuous component of the spectrum in 
billiards of class 3 was inferred by looking at scaling properties of the spectral measure, 
obtained by finer and finer numerical inversion of ([5]). In particular a multifractal 
analysis yields a nontrivial spectrum of generalized dimensions, with a Hausdorff 
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dimension Di less than 1 and a correlation dimension D 2 which rules the power-law 
decay of integrated correlations [23 [32] ■ An example is shown in figure HJb), for the 
spectral measure associated to the angle £. As we will show in following sections, the 
local scaling properties of the spectral measure are even more relevant in connection 
with transport properties. For billiards in class 3, a nontrivial scaling of the spectral 
peaks near the zero frequency is found (see figure |4jb)), in opposition to the almost 
integrable case, in which a non empty pure point spectrum is marked by not scaling 
deltalike peaks at different resolutions [24) . 

3. Higher order statistics 

3.1. Dynamical quantities and spectral analysis 

The dynamics inside the billiard table in figure [1] is equivalent to the dynamics of the 
particles in a two-dimensional infinite periodic lattice with square obstacles, i.e. in a 
generalized Lorentz gas, recently examined in [12l [13l [HI [TBI [28] . The unfolded system 
is obtained by reflecting the elementary cell and the segment of a trajectory at each 
collision point with the external square. Instead of considering particle diffusion along 
the channels of the extended system, we examine the transport generated by billiard 
trajectories revolving around the inner square obstacle. 

For this purpose, the natural observable is the angle spanned by the radius 
vector, joining the center of the billiard O with the collision point z, when the particle 
is moving from z to B v z; it is assumed positive when counterclockwise. The total 
angle accumulated by a single particle up to time t is: 

t-i 

5(z,t) = (6) 

s=0 

S(z, i) / (27r) gives the number of revolutions completed by a trajectory up to the time 
t. 

In [24[ the 2-nd order moment a 2 (t) of S(z,t), namely the r.m.s. number of 
revolutions, was examined; for generical values of the parameters in class 3, an 
"anomalous diffusion" was found, marked by an algebraic growth of a 2 in time with 
an exponent ranging between 1 (normal diffusion) and 2 (ballistic transport); this 
exponent has shown to be related to the zero-frequency scaling index of the density 
power spectrum associated to the observable £(z). In this paper we extend the analysis 
to moments of arbitrary order. 

Since £(13* z) is a "power signal", i.e. {^(B^z)} tei ^ 2 (Z), as a function of 

time t, but W — liniT^oo ^J2t=-T £ 2 {B\p z ) = C| lme (0;z) < +00, we analyse the 
trajectory of a particle on a finite time interval: —T < t < T (with T positive integer, 
i.e. T e N+). Spectral analysis will involve Fourier transform of the signal ^(B^z) on 
finite portions of trajectories. 

Hence we call £t the partial sum of the Fourier series of £(B*z): 

T-1 

i T (9-z) = (k T *g) (M = E i(K z )^ lW = Pr(P;z)eW-<'\ (7) 

t=-T 

where * denotes convolution and kx(x) = Y^t=-T e ~ lXt = e ^ (iSP^)' ^' ie u PP er 
limit of the sum in ([7]) is determined by the fact that £(B T z) is not defined. 
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Figure 3. (Colour online) Fourth order moment of E.(z,t) vs time in logarithmic 
scales, for the same data of figure[4] The slope of the straight line (= 3.4) is given 
by the theoretical prediction (13S I t with n = 4 and a(0) = 0.3. 
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is a sequence of continuous, bounded and periodic complex 

7rZ) denotes the 1-dimensional torus, 
oo of the sequence of partial sums may 



function of the frequency 8 6 T; T = 
As ^(B^z) 7^ for t ^ oo, the limit T 
not be a function in ordinary sense; as \£(t; z)\ < A Vt € Z and < A < ir, from the 

theory of Fouries series, it follows that < £t(8; z) \ converges weakly for T — > oo 

to a tempered distribution of period T. 

Since £ is a real- valued function, z) = £,t(~8; z) and the absolute value and 
the phase of £t are respectively even and odd functions of 8 (mod 2ir) : 

p T (-9;z)=p T (0;z) = p T (\e\;z) (8) 

M-0;z) = -Mbz)- ( 9 ) 

In particular £r(0;,z) = pt(0;z), 4>t{Q]z) = 0; while, from periodicity and from ([9]), 
it follows that 0t(— tt; 2) = <^t(tt; 2) = 0. 

By looking at forward trajectories starting at (s, (/?) and their time reversal, 
generated by the "backward" (inverse) operator B~ x on (s, — ip), we can verify that 
the following identity holds: 

t(B t v (s,- V )) = -aB- t - 1 (s,<p)). (10) 

This identity induces the following property of the partial sums: 

Ue;(s,-uj)) = -e i6 i* T (6;(s,uj)), (11) 

which is equivalent to pr(8; (s, —</?)) = pr(8; (s,ip)) and </>t(#; (s, — ¥>)) = + 7r — 
^r(ff; (mod (2tt)). 
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3.2. Higher order moments and correlations 

To get a more detailed picture of anomalous transport, we introduce the higher order 
statistics of S(z, t) and higher order spectral analysis of the phase variable £(z). 

While, from the numerical point of view, we may take arbitrarily real powers of 
the diffusing variable |H(«, in the theoretical analysis of sections [H H] and El we 
restrict to integer order moments of S(z, t). Integer order moments can be expressed 
in terms of multiple-time correlation functions, related to higher order spectra. 

The n-th order moment is given by: 

a {n \i)= [ dfl(z)E n (z,t) n>2. (12) 

The occurrence of anomalous diffusion, previously found for a 2 (t) [24] . is confirmed 
for moments of arbitrary orders; in figure [3] the fourth-order moment is shown, as an 
example. The crucial issue is to get an estimate for the exponents of the algebraic 
growth in time of integer moments. 

By a generalization of (j3|), we introduce a multiple-time phase-averaged 
correlation function of the observable £, which depends on (n — 1) time intervals: 

Cf(t)=/ dn(z)t(z)(l[t(B$z)), f=(f lj t 2 ,-..,t n _i)GZ B - 1 . (13) 
Jm * \i=o ) 

The n-th moment a^(t) can be expressed as a function of C^ h (i): 

t—1 t—1 t-l . 

^ (n) (*)= EE'" E / ^)^z)^z)^z)--^{B^z) 

t =0ti=0 t n ^ 1 =0^ M f 

= EE- E cf{t-t T), (u) 

t =0ti=0 t„-i=0 

where / is the vector with all entries equal to 1; note that the n-th order moment is 
related to a (n — l)-point correlation function. 



4. Higher order spectra analysis. 

4-1. Multi-dimensional Wiener- Khinchin theorem 

The (n — l)-point correlation can be expressed in terms of a (n — l)-dimensional 
inverse Fourier transform of the (n — l)-th order spectral density distribution function 
(polyspectrum) on T™" 1 [HUG!]: 

cf { t) = ^J^e m m^ (is) 

T" _1 is the (n — l)-dimensional torus, i.e. the Cartesian product of (n — 1) tori 
T x T • • • x T and 6 = (0 1 , 9 2 , • • • , 9 n - 1 ) e T"~ 1 . Equation ([T5j) is a formal generalization 
of the Wiener-Khinchin Theorem to the multi-dimensional case. However, since for 
n > 2, C^ h (t ) is not a positive definite function, according to Bochner's Theorem, 

m^(9) is not a positive distribution and m^(9)d9 cannot be interpreted as a probability 
measure on T n_1 [36j 137]. 
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An alternative version of the theorem involves the time-averaged correlation 
function: 



Cf me tfz)^lim ^ J2 Z(B» (16) 

J=l / 

d0 Si (9;z)e iff - r , (17) 



t=-T \l = l 

1 



jn-l 



with 



m c (0) = / dCl(z)s 6 {6;z). (18) 



If the system is ergodic, s^(0; z) = m^(9), for il-almost every z € M. 



4-2. Polyspectra 

We may obtain m^(9) in terms of partial sums of the Fourier series of £,(B^z). 

Firstly, we derive s^(9;z) from (|17[) . Since we consider the trajectory of a 
particle on a finite time interval, each signal £ > (B t + tl z) in ((T5)) is substituted by 
Xx T (t + ti) ■ z), where xir 1S t ne characteristic function of the interval of 

integers 1 T = {t E Z; -T < t < T - 1}. 

The square window of the signal is expressed by the Fourier transform 

XiAmB'z) = ±- f d9 i T (9;z)e M , (19) 

which is obtained by the Convolution Theorem from ([7|). 

The polyspectra s^(9; z) (and thus m^{9)) is given as a limit of a sequence of 
functions {sx.^(9; -z)}tgn + ? which converges to s^(0; z) in a weak sense: 

se(0;z) = lim stA9\z) weakly. (20) 

T— >oo ' 

By inverting (fT?]) and substituting (fT9|) . we get: 

T-l T-t-1 T-t-1 n-1 

«t,£(^) = 2tE E ••• E e- ie v T (t)e( J B^)nxi T (t+toe( J B t+t! 2 ) 

t=-T t 1 = -T-t t n -i=-T-t 1=1 

2T(^rl dd J Tni d S / iT(9:,z)k* T (0 + e(9)) f[ z)k T (9i - 00, 



i=i 



where 0(0) = SILi. 1 By using that kxix) (and /c T (x)) is an approximate identity, 
i.e. limr^oo L d9kr{9 — 9')f(9) = 27r/(0') for all continuous functions / on T, the 
final expression of s^(0; z) is obtained: 

s 6 (6; z) = T lim ^ f JJ *) j It *) • (21) 

m^(0) is derived by applying phase averages (fT%|) to both sides of (j2"Tj) . The result is 
consistent with analogous formulas in [331 13S] ■ 

We introduce the notation: Tt(9;z s j) — YYj=i Pr(Qj\ z s,i) and <&t{9',z s j) — 
S"=i <pT(0j;z s> i), with z Sii = {s,(pi) and ip t = tp - Itt/2 (I = 0,1). 
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By evaluating the phase average with ([1} and by making use of the property (fTTj) , 
we get: 

m^(9) — ^lim mx,f(^) weakly 
1 f L 1 

m T ^(9) = — ds^Jcos^ •M T)f (fi;a,, I ); (22) 



2=0 



for n odd integer: 

M T , e (0; z.,i) = 2^r T (0; ^,0/^(6(0); sin ($ T (9; z,,i) - (fr(Q(fy i ( 23 ) 
for n even integer: 

M T , e (0; = ^r T (0; z s ,Opt(6(0); z.,,) cos ($ t (0; «.,0 - <M©(#); *.,«)) • ( 24 ) 

Explicit formulas for the first few orders are given in | Appendix A| 

Higher order spectra m^{9) possess different symmetry properties. 

From (12~T1) and since ££(0; z) = £,t(—6', z), a conjugate symmetry property holds 

for m^(9): 

ml(e)=m e (-9); (25) 
this condition guarantees the reality of multiple-time correlation C^ h (i) , expressed by 

In particular, from (j23]l . (f23j> . (0) and ©, we get: 

m^{—9) — —m^(9) for odd n (26) 

m^{~9) = m^(9) for even n. (27) 

For odd n (|26l) entails that the multiple integral of m^(6) on T™ -1 is null, because 
the contributions of hyperoctants associated to 9 and —9 cancel each other. For even 
n instead, owing to (|27]l . these contributions are equal and the integration domain can 
be reduced from 2 n_1 hyperoctants to 2™~ 2 hyperoctants. 

The multiple-time correlation function C^ h (t) posseses further symmetry 

properties, which reflect to symmetries of m^{ff) and allow to further reduce 
the integration domain in (fT5|) to the "so-called" principal domain [271 132] • 
These symmetry conditions are reviewed in |Appcndix A.l| for the bispectrum and 
trispectrum, i.e. n — 3 and n = 4. 

4-3. Single frequency case: small frequency asymptotics of the spectral measure 

In section [S] we make use of the asymptotic behaviour of polyspectra for vanishing 
frequencies. Since (|2~Tj) and (1221) are expressed in terms of products of one-variable 
functions, in this section we focus on the single frequency case, namely n = 2. 

The correlation functions (|15j) and (|17|) are expressed by 1-dimensional integrals; 
s^(9; z) is the density power spectrum of the signal, obtained by a weak limit of the 
sequence of functions {st,^(9; z}t6N+ : 

*t,#; z) = ^£t(6; z) C t (0; z) = -Lp2(0; z). (28) 
By comparing (IT51) . for n = 2, with ([3]) we have 

dfi s (9) = m ( (6)d6, (29) 
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Figure 4. (Colour online) (a) Numerical reconstruction of the density power 
spectrum m^{9) associated to the angle § with a resolution A8 = 27r/2 14 , for 
= 7r/2 and tan</3 = 7r/4. (b) Mass contained in small intervals of radius 8 
centered at vs <5. The slope of the straight line gives the Holder exponent at 
= 0, i.e. a(0) = 0.30. Different symbols refer to different resolutions of m^(9), 
as explained in section [6] 



which has to be interpreted in a strict distributional sense, as in our case m^(6) will 
be a singular object. In billiards of class 3, as explained in section l2~2l the absence of 
mixing excludes occurrence of purely absolute continuous spectrum and implies that 
the correlations do not decay to zero as t — > oo; therefore, owing to Riemann-Lebesgue 
Lemma, m^{9) is not an integrable function of 9 G T. 

In [24] it was pointed out that the singularity of the measure at 9 = is essential in 
order to get an anomalous second moment of the transporting variable; in particular 
the indices that quantify such singularities are the critical Holder exponents a(9), 
which are defined for each point 9 in the support of the measure as 

nmsupMfC^)) ■ <T Q = { ° " < -!2 , (30) 
8^0 5 \ oo a > a(0) 1 ' 

where Ig(9) = [6 — 6,6 + 6] is an interval of width 26, centered in 6. The measure is 
uniformly a-H61der continuous (UaH) in an interval T&(6), centered in a particular 
point 6, if a positive constant c exists such that the mass fj,^(Ig(6)) < c6 a ^ for every 
interval l s (6) c2a(0). 

Smaller values of a correspond to stronger singularities of the spectral measure. 
The more interesting case is when < a < 1; in the following, we will consider a(6) 
varying within this range. The values a = and a > 1 correspond to a discrete and 
absolutely continuous component of the spectral measure, respectively; if a(6) = 0, 
A*?({f?}) > and if a(9) > 1, the measure is continuous and differentiable in Xa(6) 
and the derivative m^(9), being an integrable function on T&(9), is the density of the 
measure, with respect to Lebesgue measure. According to numerical approximations 
of the spectral measure, billiards belonging to class 3 are characterized by values of a 
in the range < a < 1; in figure |4] a typical case is shown, in which a(0) = 0.3. This 
is consistent with the presence of a singular continuous component of the spectrum; 



Higher order statistics in the annulus square billiard 



13 



however, an exponent a € (0, 1) is not a sufficient condition to have either a continuous 
[24] or singular continuous part 25 of the measure. 

In section [5] we will use a relationship which connects the Holder exponent of the 
measure at some point 9 to the asymptotic behaviour as T — > +oo of the sequence 
{wit,£(0)}tgn + ; it is based on the equality m,T,$(0) = ^ Jjdfi^(9')KT(0 — 9') with 
Kt{x) = (smTcc/sin(x/2)) 2 /2T. This relationship is derived in [33] and reviewed in 
|Appcndix B| 

If the measure is UaH in an interval T&(6), then there exist a positive constant 
D and T = T(A, 9) such that for T > f : 

m T ^{9) < DT 1 -^ uniformly for 9 G Z&{9). (31) 

As we are dealing with an ergodic system, the bound (|3"Tj) holds also for st,^(9; z), 
for fi-a.e. z because the dependence on z in (jT5J) is actually missingU 

Consequently, it follows from (|28l) that, for f2-a.e. z, the sequence of functions 
pr{9; z) as T — » +oo is bounded by: 

p T (9; z) < CT 1 "^ for 9 e X A (0), with C > 0; (32) 

in particular, for instance, if 9 = and A = 2ir/T (with T > 2), (|32| holds in every 
interval I 5 (0) with <5 < 2tt/T. 

5. Analytic estimate for the moments' scaling function 

An accurate study of moments' asymptotic behaviour at large times is the central 
point of the present paper. The moments' scaling function j(n) is defined as the real 
number j(n) such that the discrete Mellin transform: 

+oo 

/(n) ^) = E TIT? ( 33 ) 
t=l 

converges for (3 > 7(71) and diverges for < 7(n). This definition of 7(n), respect 
to other definitions based on the asymptotic behaviour of a^ n '(t) as t — > 00, has the 
advantage that it does not take into account of possible subdominant contributions to 
the transport process. 

The moments' scaling function may be written as -f(n) = v{n)n. Normal 
transport and weak anomalous transport fall in the category with v(n) = vq = const. 
with Vq = 1/2 and v$ ^ 1/2, respectively; at long times, the distribution of S(z,i) 
relaxes to a self-similar function, which is a Gaussian distribution when = 1/2. 
Strong anomalous transport corresponds to the case where the distribution of S docs 
not collapse to a self-similar form, multiple scales exist and a phase transition is 
observed, marked by a piecewise linear function v(n) [29l [33] . 

For even moments of integer order ( n > 2), we obtain an analytic relationship, 
linking the scaling function 7(71) with Holder exponent at 9 — of the spectral measure 
associated to t;(z). This relation extends to higher order moments the formula found 
in [21] for the exponent 7(2) of the 2-nd moment. 

X More precisely, accordingly to 1251 . the critical exponent 1 — a(9) in (13 1 1 is an upper bound for the 
exponent of sx^(8', 2), for f2-a.e. z £ Ai v . 
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By substituting (fT5|) into ([TO , the moment of order n is expressed by the following 
multiple integral over the (n — l)-dimensional torus T™ : 

^ (*) = TTl^T / d0 m s (0 )A (0(0)) J] <^ ) 

where D t (x) denotes the kernel: 



= e-l^*" 1 ) ^ e — = (35) 

n=0 Sm ( 2 , 



(f*) 



The conjungate symmetry property ([2"S")) of the polyspectrum guarantees that the 
expression (j34j) is real. Moreover, since D t (x) is an even function of x, the parity 
transformation property of the integrand in (|34[) is determined by m^(9). From (|26p , 
it follows that all moments of odd order are null. This is consistent with the fact that 
the angle H(z,i), accumulated by a single particle, may assume positive or negative 
values; hence, in |Appendix C| odd order moments of |S(z,t)| are taken into account. 

For even n, owing to symmetry (|27[) . the integration domain of (|34l) can be reduced 
to half of the (n — l)-dimensional hyperoctans. 

By substituting (l34l) into (j33l) . we get: 

J, (n) (/3) = / d6m li {d)S{p,e)= lim / d (9 m r ,£ (<?)5r(/3, 0) (36) 

r 1 n-l 
5 T (£, (9) = £ ™ A (0(0)) J] 508, 0~) = ,Hm 5 T (/3, 9) 

t=i j=i 

The integration domains Lj 1 in (|3^|) are Cartesian products of (n — 1) half tori, 
T + = [0,7r] or T_ =] — 7r, 0]. For instance for n = 4, are the following octants: 
Lf = T+ x T+ x T+, L| = T+ x T+ x T_, L| = T+ x T_ x T+ and Lf = T_ x T+ x T+. 
As shown below, the convergence of ((3) is guaranteed under the condition 

0>n(l-ffl); (37) 

a(0) is the Holder exponent at of the spectral measure, associated to £. 

Therefore, according to the definition of the moments' scaling function j(n) 
defined via the discrete Mellin transform (|3"3"|) . we have: 

7(n) = n I 1 — I n positive even integer. (38) 



As stated by (f38|) . the spectrum of moments is governed by a single scale, and no 
strong anomalous transport takes place. 



5.1. Derivation 0} the convergence condition |37[ ) 

The argument (for even order moments) consists of two different steps. The first 
step involves the second term appearing in ([55]). namely St{/3,0), which constrains 
an evaluation of (|36p around the origin. Indeed, lim^oo D t (x) = uniformly 
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for | a; | > 8 > and the dominant contribution to the integral (I36[) comes from 
simultaneous zeros of the denominators of D t (9j), namely from Oj = 0, Vj = 1, • • • n— 1 

M- 

The kernel D t (x) is an even function of x, whose smallest positive zero is given 
by x = 2ir/t; moreover, in the interval |x| < 27r/i: < D t (x) < t. It is sufficient to 
choose \9\ < 27r/((ri — l)T) to get all the arguments of the kernels in St((3, 9) within 
the range |a;| < 2ir/t. 

Around the origin, we may thus write 

T 

ST(/3,fl1<£^i"^(l,T), for 0<_??L_ (39) 

For P <M, the partial sum is bounded by: 

s(l,T) < c{fi)T n -^ (40) 

where c(/3) is finite function for j3 ^ n. 

By using T < 2ir/{{n - 1)\6\), the final result for S((3, 9) is: 

S(p,6)<s(l,T)<ci(p)\9f- n (41) 

where: ci(/3) = (2vr/(n - l)) n -Pc(0). 

Secondly, we derive the local scaling behaviour of m^{9) near a singularity 9 
in frequency space T Tl_1 . As in 1-dimensional case, the scaling is related to the 
asymptotic growth of the sequence {mx^ (0)}t6N + as T — > oo. This link holds locally 
in frequency space, inside a small (n — l)-sphere centered at the singularity of radius 

\S9\ = \9- 9\ < 2ir/(T(n- 1)); then for large times, i.e. T -> oo, \59\ -> 0. 

The high-dimensional spectrum (|2"2")l is expressed by a phase average of the 
contributions of different trajectories, i.e. Mt,£ (9; z s ,i)- Formula gives Mt^ as a 
function of px and </>t of a single variable 0, where 9 means 9j with j = 1, • • • , n — 1 
or 0. By making use of the asymptotic behaviour (|32[) of the modulus px, valid for 
T — > oo, from we get the bound: 

n-i 

|M r , c (ei, • • • , e„-i; < C"T"- 1 T- a (®)/ 2 J] T" a ^)/ 2 , (42) 

i=i 

and, in particular, in a small (n — l)-sphere centered at 0, 

m^e) < C\9\<^- l ) +l for \9\ < j^f- (43) 

In lhs of (|43p we omit the absolute value because, as explained in section [31 for a fixed 
T, the modulus pt and the phase 4>t are continuous functions of 9 and in particular, 
since 4>t(9 = 0) = 0, 4>t — > for 9 — > 0. Hence the scaling behaviour in of the total 
phase of Mt,£ , which is a sum of the single phases, is a trivial. 

We finally derive condition (|37|). under which the integrals (f36| converge. 

We consider a (n — l)-sphere of radius 2tt/((ti — 1)T), centered in the origin of 
the frequency space, inside which the integrand function is bounded by (|4T]) and (f43j) . 
Dominant contributions to integrals (I36|) are restricted inside this (n— l)-dimensional 
sphere. 

The integrals may be evaluated by (n— l)-dimensional hyperpherical coordinates 
{r, ^1,^2, • • • 01]. The vector 6* may be written as 9 = rw; < 



§ For p > n the integral J 36 1) is always convergent in a small sphere around the origin. 
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Figure 5. (Colour online) Scaling function 7(9) of the moments cr^ q '(t), given 
by (145 b . plotted vs the order q. Dotted straight lines have a slope q ^1 — ° 2° y ' 
according to the theoretical prediction H38H , extended to real values of q. Moments 
of integer orders are marked by halos. (a) refers to same data of previous figures, 
for which the scaling exponent o(0) = 0.3 is derived in figure [4] The parameters 
of (b) arc: h/l 2 = l)/2, tan ip = it/ 4 and a(0) = 0.47. 



r < 2n/((n — 1)T) and Q is a unity vector with components: u)\ = 
cosipi,uj2 — sin -0! cos ^2 ; w ri-2 = sin 4>i sin ip2 • • • sin 4>n-3 cos VV1-2 ; = 

sin^! sin -02 ' ' 1 sin-0„_ 2 with < 0/ < 7r for / = l,n — 3 and < ip n -2 < 27r. The 
Jacobian of the coordinate transformation is J = r™ -2 sin™ -3 ipi sin™ -4 ip2'" sin-0 n _3. 

The integration over variables r and uj can be carried on independently. From 
(|41l) and (|43|) . we may write: 



4 n) (/3) „ 2 <c 2 (/3)/ duj drr' 9- "( 1- ^ 1 ) -1 (44) 

s n _2 is a (n— 2)-dimensional hypersphere of unitary radius and e 2 (/3) is a finite function 
for (3 ^ n. 

The integral over the angles is trivial, while, the radial part is convergent under the 
condition (l37l) : /3 > n(l — a(0)/2), which finally yields the estimate for the exponent 
of even order moments. 



6. Numerical simulations and conclusions 



From the numerical point of view, we may extend the analysis to moments of the 
diffusive variable |H(z,t)| of arbitrary positive real order. The q-th order moment is 

a^(t) = [ dn(z)\E(z,t)\ q , qeR+; (45) 
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the scaling function j(q) is defined by (|3"3")l . with <jW replaced by a^ q \ Note that in 
spite of having introduced the absolute value respect to the definition (|T^|) . we keep 
the same notation. 

Numerical simulations, presented throughout the paper, refer to a billiard table 
belonging to class 3 with parameters lx/l% — tt/2 and v x /v y = tan tp = 7r/4; a second 
case is shown in figure E](b), referring to h/h = (%/5 + l)/2. 

The analytical expression for the exponent 7(71) of the algebraic growth of integer 
order moments, given by (|38j) . is numerically tested in figure El in which the A-th 
order moment is plotted as a function of time t in logarithmic scales. The straight line 
has a slope 7(4) = 3.4, given by the formula with n — 4 and a(0) = 0.3. Numerical 
procedure to derive a(0) is explained below. 

Figure [5] displays the scaling function 7(9) of absolute moments of arbitrary real 
order (|43| . for the two parameters' pairs (h/h, tan f) specified above. Numerically 
evaluated exponents are shown by dot symbols as a function of q; they are calculated 
by a linear least-square fit of \na^ q \t) versus Int. 

Analytical estimates of the exponents of integer order moments are shown by 
halos, surrounding the dots. Inside the inspected range, numerical data are consistent 
with a linear dependence of 7(g) on q. The ratio r y(q)/q appears to be constant for fixed 
parameters (h/h, tan^) and to depend only on the spectral index a in 9 — 0. This 
ratio ranges between 1/2 (normal diffusion) and 1 (ballistic dispersion). Therefore 
the square annulus billiard displays a so called "weak" anomalous diffusive (or strong 
self-similar) process [29| [30] . 

We finally provide a few further details on the numerical procedures. 

The discrete-time directional dynamics at fixed ip is evolved up to a time T = 2 20 
and the number of initial conditions employed in phase averages ranges from 10 to 
10 6 . 

An approximation of the microcanonical measure inside the phase space M.^ is 
obtained as follows. The phase average is evaluated by taking a uniform distribution of 
the particles inside the accessible region of the billiard and the sign of the components 
of the unity velocity vector v is assigned at random. The first collision points with 
the boundaries of the billiard are taken as initial conditions of the Birkhoff mapping. 

Spectral analysis of the signals is reconstructed numerically by employing the 
Fast Fourier Transform algorithm (FFT), i.e. a finite discrete Fourier Transform for 
vectors in C 2T . Two different methods have been tested to calculate the coarse- 
grained approximation of the (1-dimcnsional) average density power spectrum (|18p . 
According to (fT5)). m^(9) may be approximated by a direct-FFT of the phase-averaged 
correlation sequence C^ h (t). Alternatively, according (fT8|) and (|28|). we calculate the 
square modulus of the direct-FFT of the signal ^(B t z) of single trajectories and then 
make a phase average over the initial conditions. Owing to the finite time interval 
— T < t < T, we multiply time sequences, C^ h (t) or S s (B t z), by a proper windowing 
function |42j ; to reduce the amount of negative values in the former method, we prefer 
to use a triangular window function Wt = 1 — instead of a square windowing 

function as in the theoretical treatment, because its partial Fourier series is a non 
negative function. In both methods, the resolution in frequency space is A.9 = tt/T. 
For the finest resolution, i.e. T — 2 20 , the order of magnitude of negative values, in 
the first method, is < 10~ 4 . Apart from numerical errors, the two methods give the 
same results, for a fixed resolution. By comparison of the two methods, we get that 
less than the 1% of the total mass on the torus is affected by numerical errors. In 
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figure0](a) a reconstruction oim^(9) is shown, calculated with a resolution of 7r/(2 14 ). 

The first two generalized dimensions Di (information dimension) and D 2 
(correlation dimension) [43] are calculated by making a sequence of dyadic partitions 
of the coarse-grained approximation of m^(9). At step N, the interval (— 7r,7r] is 
divided into 2 N sub-intervals Ijyj (j = 1,---,2 N ) of length 5n = 2n/2 N , each of 
which contains a mass /i(Jjvj). As N —> oo (i.e. 5n — > 0), D\ and D 2 are defined by: 

Xi(N) = J~] i" (I N,j ) In M at, j ) ~ -Di ln5 w 
3=1 

X2 (A) = ln^/i 2 ^) - #2 In (5a,. 

i=i 

We extrapolate D\ and D2 by a linear fitting in logarithmic scale with N ranging from 
N = 2 to N = 20; an example is shown in figure [5] (b). A value of D x < 1 is found, 
consistent with the presence of a singular continuous component of the spectrum. As 
shown in figure [5] (a) , the obtained value of D 2 gives a satisfactory estimate of the 
power law decay exponent of the integrated correlation function. 

The scaling index a of the spectral measure in 9 = is also evaluated numerically. 
In figure HI (b) the mass contained in the interval [— <5, S] is shown as a function of 5. 
Intervals of different widths 5m = M A9 are considered, with a resolution A9 ranging 
from 7r/2 20 to 7r/2 6 ; different symbols refer to M — (triangles), M — 1 (squares) 
and M — 8 (circles). The mass is given by 6m, 5m]) = Ylj=-M m iU^) anc ^ ^ 
is supposed to scales as ~ <5"/°' ) . The straight line, which fits the data in logarithmic 
scales in the range —11 < ln6 < —7, has a slope of a(0) = 0.30 ±0.05. The data in the 
left part of figure H] (b) are not taken into account in the derivation of a(0) because 
the corresponding mass is close or below the numerical precision of the simulations. 
The same calculation in case h/h = l)/2 gives the value a(0) = 0.47 ± 0.05. 

Note that the error on the numerical estimate of 5(0) is not relevant when compared 
with data of figure O 

We have thus provided analytic arguments and numerical support that polygonal 
billiards may exhibit anomalous transport in a weak sense |29j : the distribution of 
the angle, accumulated by single particles at fixed time, is not gaussian and the 
moments' asymptotic growth is governed by a single scale, which is related to the 
Holder exponent of the spectral measure at zero. 

This work has been partially supported by MIUR-PRIN project "Nonlinearity 
and disorder in classical and quantum transport processes" . 

Appendices 

Appendix A. Power spectrum, bispectrum and trispectrum 

For n = 2, namely for a single frequency 9 € T, the result ([28} is retrieved from (|24|) : 
the weak limit s^(9; z) of the sequence (|28|) is the density power spectrum of the signal. 

For the first few orders, i.e. n = 3 and 4 (bispectrum and trispectrum), we have 
the following explicit expressions, respectively: 

M T ,z(9i,9 2 ;z s .i) = —pt(&i; z s .i)p T (9 2 ; z s .j)p T (9 1 +9 2 ;z s j) ■ 

■ sin(£r(0i;«,,j) + <t> T {0 2 ; z s ,i) - <£r(0i + h;z s ,{)), (A.l) 
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and 

M T ^(d x ,0 2 ,9 3 ;z s j) = ^=Pt{Qi\ z s j)p T (6 2 ; z s> i)p T (9 3 ; z s .i)p T (6 x + 2 + 6 3 ;z s j) ■ 

• cos(0 T (6>i; z s ,i) + <p T (0 2 ; z s ,i) + <p T (0 3 ; z s ,i) - fa {61 + 2 + 3 ; z s> i)), (A.2) 
with z s j = (s, ipi) and (pi = ip — lnj2 (I = 0, 1). 

Appendix A.l. Symmetry properties 

Owing to the symmetry properties of higher order spectra and of multiple-time phase 
averaged correlation functions (|13l) . the integration domain in (|15l) can be reduced. 
Some properties, as (g5|), ([HI) and ([27]), follow directly from (gT]) and ©, © and have 
been used to restrict the domain in (|36[) , as explained in section [5] 

A general procedure to determine the principal domain of polyspectra, i.e. the 
nonredundant region of computation, is explained in [351 (Ml El]- The symmetry 
properties of higher order spectra can be grouped into: 

(i) invariance under permutation of any pair of frequencies; 

(ii) the conjugate symmetry property (|25p. which implies the redundancy of half of 
some frequency axes; 

(iii) periodicity of period 2tt is each frequency variable; this is a consequence of the 
fact that the partial sums (J7J are periodic functions of the frequencies 61 € T 
(i = l,---,n). 

In this appendix, we review some symmetry conditions for the 2- and 3- dimensional 
cases [27] . 

Two dimensional case. From the definition of 2-time phase-averaged correlation 
function, i.e. ([TBf with n = 3, and the invariance of the measure dtt(z) under B^, we 
get: 

C ph (t x ,t 2 ) = Cv h {t 2 M) = C^i-hM - t x ) = CP h (-t 2 ,t x - t 2 ). (A.3) 
From (jA.ll) . © and ©, we have: 

m(d x ,0 2 ) = m(6 2 ,&x) = m*(-0 x , -6 2 ) = m(6 x , -6 X - 6 2 ) = m(9 2 , -6 X - Q 2 ). (A.4) 

Owing to (|A.4|) . the bispectrum is symmetric about the lines 9 X = 9 2 , 0\ + 29 2 = 2im x 
and 9 2 + 29 x — 2ixn 2 (n x ,n 2 € Z). The principal domain is the triangle of 
vertices A(0,0), B(ir,0) and C(§tt, §tt), i.e. < X < |tt, < 2 < 6 X and 
§7r < 9 X < 7T, < 6 2 < -29 x + 2tt. 

Three dimensional case. The 3-time phase-averaged correlation function satisfies: 
C ph (t x ,t 2 ,t 3 ) = C ph (t x ,t 3 ,t 2 ) = C ph (-t 1 ,t 2 ~ t x ,t 3 - t x ); (A.5) 
and the trispectrum: 

m(d 1 ,9 2 ,6 3 ) = m(0 x , d 3 , 6 2 ) = m(-9 x , -d 2 , -6 3 ) = m(9 2 ,9 3l -9 X - 2 - 6 3 ). (A.6) 

Cyclic relations, obtained by permutations of (t x , t 2 ,t 3 ) in (IA.5I) and (9 x ,9 2l 9 3 ) in 
(IA.6|) . hold. The trispectrum is symmetric about the plane 28 x + 6 2 + 6 3 = 2-Kn x (and 
cyclic). The principal domain is a polyhedron with vertices A(0, 0,0), £?(?, 77, 5), 
COr, 0,0), D(§tt, §tt,0), B(7T,0,-f) and F(7r,7r,-7r) [3J]. 
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In this appendix we sketch how to derive (|31[) . 

We assume that the spectral measure dfi^(9), associated to £, is uniformly a- 
Holder continuous (UaH) in an interval 1a(9) = [9 — A, 9 + A], namely a positive 
constant c exists such that 

j.0+6 _ _ 

^{l s {9))= d^(9')<c8 a{e \ Vls(6)cX A {6). (B.l) 

J6-5 

The sequence {tot,£ (^)}tsn + is defined by (fTg)) and (|2"51) as 



m T ^(9) 



<m{z)s TX {6; z) 



(B.2) 



= 7^ dn(z)t T (6;z)&(6;z) 



By making use of ([7]) and then of |3]), (TT5l) and (|29|) , we obtain 

T-1 T-1 : 

' d^(9')K T ( 



i=-T t' 



1 

2^ 



(B.3) 



with 

*k«) = ^ E E 



T-1 T-1 



2T 

1 /sin To; 



ii=-Tt 2 =-T 
2 



2T 



(d T (x) - 2d T (a;) cosTz + l) 



2T V sin f 



(B.4) 



Kt(x) is the Fejer kernel, up to a constant factor, and dx{x) is the Dirichlet Kernel: 

(B.5) 



t=-T 2 



The kernel (|B.4|) is a positive, periodic even function of x € T; moreover, it fulfills: 
(t) K T {0) = 2T; 

ii) lim Kt{x) = uniformly Ixl > 5 > 0; 

T— >oo 



m) < ifr(x) < 2T \x\ < 



T 



T 
2~k 2 



< \x\ < (k + 1)^ and k > 1. 



iv) < K T (z) < 2// > 

«) ^<^t(^)<2T H<^. 

By using the properties above, the proof of theorem 2.1 in |25] can be reproduced 
for rriT.^(9)- Note that rriT^{9) corresponds to Gt(9), and (3 to 1 — a in |25] , 

Let £ Ia(9) and < <5 < 7r, we show that (pTTj) holds in an arbitrary 
interval Is(9) C 1^(9). X-s(9) can be covered by a finite sequence of intervals: (1) 
2T(°) = {0' G T/ |0 - 0'| < f }, in which (m) holds, and (2) = {9' G T/ fef < 
|0 - 0'| < (fc + l)f < in which (iv) holds. 
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Moreover, assumption (IB.ip implies that the mass in every interval of length n/T, 
included in T&{6), is bounded by ^{1^_(9)) < DT~ &{ -®\ with constant D > and 
T > f (A, 9); T is chosen so that for T > f : (9) C 1 A {0). 

Therefore, 



mr,dO) \,m - I {/ i(o) +E/ IW }^(^(^ - 

/ [5T/7T] \ 

< C*i + C* 2 ^ p r -a(0) T < (B.6) 



k 2 
fe=i 



denotes the integer part. 



Appendix C. Upper bound for absolute moments of odd order 

As mentioned in section [5l the moments of odd order of the observable S(z,t) vanish 
due to the symmetry property (|26p of m^(9); we therefore consider integer moments 
of the absolute value of the observable t): 

*W(t)= [ <m(*)|H(*,t)l n . (ci) 

For even n (|C.1[) reduces to (fT2")) : we are interested to n odd integer (n > 1). 



Since |H(z,i)| < YllJo \£(B^', Z )\, * ne following bound holds: 

<*<">(*)= / dO(z)fV|e(B s z)|^ . (C.2) 



Therefore, the exponent 7(n) of the algebraic growth of <j( n )(i) constitutes an upper 
bound for the exponent 7(71) of (|C.1|) : 7(11) is again defined via a discrete Mellin 
transform (|33[) . with er replaced by cf. 

We reproduce calculations in section [4j by starting from the new phase 
observable ^(B t z) = \£(B f z)\. The partial sum of Fourier series, i.e. £ T (9;z) = 
Y%=-T ^{B* z)\e- ie \ verifies the property (cf. (HI])): 

t T (6; (s, -?)) = +e w T T (9; (s, <p)). (C.3) 

By using (IC3|) in the phase average, the following expression for Mt^ is found (cf. 
(EH and (HQ)): 

M T ,M = T (9; z s ,i)pr(e(9); z s ,i) cos ($ T (0; - 6r(9$; «.,,)) . (C.4) 



Mt,£ is now an even function of 9 and it shares the same symmetry properties of (|24|) . 
For n = 3, in particular we have: 

Mt^(9i,9 2 ;z s j) = — Pt(8i; z s ,i)pt(02; Zs,i)Pt(9i + 9 2 ;z s ^) ■ 

• cos(0r(ffi; z„,i) + $ T {0 2 ; z„,i) - 4> T (0i + 2 ;z s> i)). (C.5) 

The exponent 7(71) can be derived by proceeding as in section[5j by denoting a(0) 
the Holder index of d/iifi in 0, we get 7(71) = n(l — S(0)/2). 
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The mass associated to |£| in a small interval, centered in 0, is 

T — 1 T— 1 

M|||(X 5 (0))= lim A f dfi(s) £ ^ |£(B>)£(B^)|sinc(5(i-i')), (C.6) 

Jjvi, p t=-T t' = -T 

with sinc(a;) = sinx/x. If S < n/(2T), ^(ls(0)) < ^i(ls(0)) and therefore 
a(0)<a(0). 
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